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Abstract 

We review recent developments regarding non-equilibrium quantum dynamics and many-body physics with 
light, in superconducting circuits and Josephson analogues. We start with quantum impurity models addressing 
dissipative and driven systems. Both theorists and experimentalists are making efforts towards the characterization 
of these non-equilibrium quantum systems. We show how Josephson junction systems can implement the 
equivalent of the Kondo effect with microwave photons. The Kondo effect can be characterized by a renormalized 
light-frequency and a peak in the Rayleigh elastic transmission of a photon. We also address the physics of hybrid 
systems comprising mesoscopic quantum dot devices coupled to an electromagnetic resonator. Then, we discuss 
extensions to Quantum Electrodynamics (QED) Networks allowing to engineer the Jaynes-Cummings lattice and 
Rabi lattice models through the presence of superconducting qubits in the cavities. This opens the door to novel 
many-body physics with light out of equilibrium, in relation with the Mott-superfluid transition observed with 
ultra-cold atoms in optical lattices. Then, we summarize recent theoretical predictions for realizing topological 
phases with light. Synthetic gauge fields and spin-orbit couplings have been successfully implemented with 
ultra-cold atoms in optical lattices — using time-dependent Eloquet perturbations periodic in time, for example — 
as well as in photonic lattice systems. Einally, we discuss the Josephson effect related to Bose-Hubbard models 
in ladder and two-dimensional geometries. The Bose-Hubbard model is related to the Jaynes-Cummings lattice 
model in the large detuning limit between light and matter (the superconducting qubits). In the presence of 
synthetic gauge fields, we show that Meissner currents subsist in an insulating Mott phase. 


Keywords: Condensed-Matter Physics with Light, Superconducting Circuit Quantum Electrodynamics Networks, 
Josephson Effect and NanoScience, Dissipative and Driven quantum impurity Models, Jaynes-Cummings and Rabi 
lattices, Topological Phases and Synthetic Gauge Eields, Eloquet theory 


1. Introduction 

The field of quantum electrodynamics [1] has attracted some attention with novel and modern applications in 
non-equilibrium quantum optics and quantum information [2]. In fact, many experiments have accomplished to 
couple light and matter in a controlled manner using cavity quantum electrodynamics. Quantum matter here can 
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be either Rydberg atoms (cold atoms) or trapped ions [3, 4] for example. A step towards the realization of many- 
body physics has also been made through the realization of model Hamiltonians such as the Dicke Hamiltonian [5], 
where the associated super-radiant quantum phase transition has been observed in non-equilibrium conditions [6]. 
A solid-state version of cavity quantum electrodynamics, related to circuit quantum electrodynamics, built with 
superconducting quantum circuits [7, 8], is also a very active field both from experimental and theoretical points 
of view. Theorists have predicted novel emergent quantum phenomena either in relation with strong light-matter 
coupling [9] or non-equilibrium quantum physics [10, 11]. 

The goal here is to review developments, both theoretical and experimental, towards realizing many-body 
physics and quantum simulation in circuit QED starting from small networks to larger ensembles of superconduc¬ 
ting elements in the microwave limit. An experimental endeavor has been accomplished towards the realization 
of larger arrays in circuit QED [12, 13, 14, 15] and towards controlling trajectories in small systems [16, 17]. 
This research is complementary to existing efforts on collective effects in ultra-cold atoms [18, 19, 20, 21, 22], 
Nitrogen-Vacancy centers [23] and opto-mechanics for example [24, 25]. We shall also summarize (our) progress 
related to quantum impurity physics with light, such as Kondo physics [26, 27, 28, 29, 30, 31] of photons 
[32, 33, 34, 35] and non-equilibrium quantum dynamics [36, 37, 38, 39, 40], as well as theoretical predictions 
concerning many-body physics in larger network circuits described by the Jaynes-Cummings or Rabi lattice 
[41, 42, 43, 44, 45, 46, 47, 48, 49]. Some progress has been done towards the implementation of an effective 
chemical potential for photons [50]. The effect of dissipation will be discussed for small and larger networks. In 
Ref. [51], we have discussed other theoretical aspects in relation with these quantm impurity models; see also Refs. 
[52, 53]. It is also important to mention recent theoretical developments regarding the strong-coupling light-matter 
coupling, both at equilibrium [9, 48, 54, 55, 56, 57, 58, 59] and in non-equilibrium [36, 60]. Artificial gauge 
fields for photons can be realized through nano-circulators, consisting of superconducting Josephson-junction 
rings penetrated by magnetic fluxes for example [61, 62, 63], allowing the realization of topological phases [64]. 
The stability towards disorder effects was also studied [64]. An experimental effort to build such nano-circulators 
has been realized in Ref. [65]. In fact, inspired by the discovery of quantum Hall physics [66, 67, 68] and 
topological insulators [69, 70, 71], experimentalists have already succeeded in realizing topological phases of 
light and spin-orbit type coupling [72, 73, 74, 75, 76] (but not yet in microwave superconducting circuits), for 
example through Floquet time-periodic drive (for recent developments, see Refs. [77, 78]); some aspects of 
topological phases with light have been discussed in the recent reviews [79, 80]. It is also important to underline 
the experimental progress in the field of ultra-cold atoms to realize synthetic gauge fields [81, 82, 83, 84, 85, 86]. 
Haldane and Raghu made a pioneering contribution in this field [87] that was soon confirmed experimentally 
[72]; they have shown that one can realize quantum Hall phases of light through magneto-optical coupling or 
Faraday effect. The (chiral) transport of photons at the edges in these systems has been measured [79] by analogy 
to quantum Hall systems [88]. This opens the door to novel applications in quantum information and quantum 
computing by analogy to electron systems [89]. The direct measurement of topological invariants [90] in photon 
lattices is in principle also feasible [64, 91, 92, 93] in particular, using the concept of anomalous quantum velocity 
[94] which has been developed further for photon systems [95]. A recent work by the MIT group on photonic 
gyromagnetic crystals with ferrite rods reports the measurement of a Chern number with photon systems [96]. 
Some aspects of the Berry phase [97] and Chern number [98] are accessible in small superconducting circuits 
[99, 100, 101, 102]. Artificial graphene with Dirac photons [103] and three-dimensional Weyl systems [104], as 
well as analogues of flat band models [105], have been realized recently. 

Other photonic arrays, including quasi-periodic arrays [106] and one-dimensional geometries [107, 108], 
have been envisioned and some have been realized experimentally [109]. 


2. Small Circuits 

This Section is organized as follows. First, we will introduce the physics of small superconducting circuits. 
Then, we will describe progress on quantum impurity systems, at and out of equilibrium. We will also account for 
dissipation effects and show that it can produce novel many-body physics, for example, based on vacuum quantum 
fluctuations associated with long transmission lines. We will describe the Kondo effect of light [32, 33] and 
address recent developments in hybrid systems [110], combining circuit QED and mesoscopic electron quantum 
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dot systems [111, 112, 113, 114, 115]. Note that in relation with quantum vacuum fluctuations, the dynamical 
Casimir effect was observed in superconducting circuits [116]. 

2.1. Superconducting Elements 

These superconducting circuits embody simple building blocks for quantum simulation, based on harmonic 
oscillators and artificial atoms - here, two-level mesoscopic systems - such as charge, flux or phase qubits. 
Nowadays, qubits in superconducting circuits yield a long decoherence time that can reach a few micro-seconds. 
For a recent progress on superconducting qubits, see for example Refs. [117, 118, 119, 120, 121, 122, 123, 124]. 
Microwave photons can propagate in on-chip resonators which can be identified to Fabry-Perot type cavities made 
from finite sections of transmission lines. The physics of such transmission lines is reminiscent of an LC (harmonic) 
oscillator described by the raising and lowering operator a'' and a. The single-mode approximation is usually 
appropriate for such cavities of t)q5ical size ~ 1cm (operating in the microwave domain) [125, 126]. At a simple 
level, (a-Fa^) represents the voltage fluctuations in the cavity whereas {a — a'')/i mimics the conjugate flux variable 
[125, 127]. Superconducting materials used in such devices are generally Niobium or Aluminum. The propagation 
of micro-wave photons in the system will be described below by coupling the cavity to transmission lines (see 
Sec. 2.3). The input AC signal Vq- cos(a)t) will be related to the time-averaged input power = Vq/[2R) where 
R is the resistance of the (left) transmission line. This will produce a mean number of photons N such that 
Pin = Nhco and N = dN/dt. The light-matter coupling will produce non-linearities in the harmonic oscillator 
spectrum. The superconducting circuits allow a great flexibility afforded by the nature of the nano-fabrication 
system; parameters in the Hamiltonian can be widely tuned with conventional lithography. Qubit frequencies can 
be set from typically 2 to 15 GHz and couplings can range from a few KiloHertz to nearly IGHz. As the lateral 
dimension of the transmission lines is typically a few micrometers (which is much smaller than the wave-length), 
the electric held in the cavity is relatively large and a strong dipolar coupling can be achieved. Gapacitors at the 
ends of the cavity control the photon leakage rates and in multi-cavity systems, the hopping rates between nearest 
neighbor cavities [10]. 

A circuit containing the cavity, artificial qubit, dissipation and AC drive effects is drawn in Fig. 1. 

2.2. Rabi, Dicke and Jaynes-Cummings Models 

The model describing the strong-coupling between a qubit and the electro-magnetic held of the cavity 
corresponds to the Rabi model (1936) [128] : 



( 1 ) 


Here, coq represents the frequency of the quantized harmonic oscillator (quantized electro-magnetic held), A 
is the resonant frequency between the two levels forming the qubit and CT; where i = x,y,z are the Pauli spin 
operators. Furthermore, g represents the dipole coupling between qubit and photon. If represents the two 
charge states of a Cooper pair box [129, 130, 131] then the capacitive coupling between a qubit and a resonator 
realizes a quantum impurity model of Rabi type [132]. An analogous model is obtained in the case of an inductive 
coupling between the qubit and the cavity. 

Although the Rabi model represents one of the simplest interacting quantum models, it does not correspond to 
a simple theoretical problem. Nevertheless, recently analytical solutions have been pushed forward in connection 
with the experimental possibilities in superconducting circuits, and more generally in nano-circuits. In particular, 
the limit g/too ~ 0.1 is commonly reached experimentally [133, 134]. The discrete Z 2 symmetry associated with 
the Hamiltonian (a ^ —a and cr^ —o^ (cj^ ^ —ct^)) cannot be ignored in this regime making analytical 
solutions more difficult to And. In this respect, it is important to emphasize the theoretical progress first realized 
by D. Braak in the strong-coupling limit [9, 54, 55, 56, 57, 58, 48, 59], concerning the solvability of the Rabi 
model. It is also relevant to underline other theoretical progress on the integrability of related quantum impurity 
models with light [135, 136, 137]. One relatively ‘simple’ limit of the Rabi model in the strong-coupling limit 
is the so-called adiabatic limit [138, 139], corresponding to the highly detuned limit A ^ 0, which has been 
addressed for example in Ref. [57]. In Refs. [37, 38], in collaboration with R E Orth and A. Imambekov, one of 
the authors has developed a relatively simple stochastic scheme to address certain non-equilibrium properties 


3 




Figure 1 : (Color online! Example of Superconducting Circuit Quantum Electrodynamics, with an artificial atom with two levels (spin-1/2 or 
qubit). The electromagnetic resonator (cavity ~ 1cm) is described through a one-mode transmission line and is capactively coupled to an AC 
drive (input signal). We show an example of spin dynamics in the strong-coupling light-matter regime. We apply the methodology described in 
Ref [36] and consider the far detuned or adiabatic limit of the Rabi model: g/too from 0.7 (green) to 1.5 (yellow) and A/coo = 0.05. The 
system is prepared in the initial state | -I- z) (gi (0) = |0, -l-z) and the revival probability of the the initial state can be computed exactly in this 
limit [140]. The positions of the maxima only depend on the cavity frequency coq- We also show the Jaynes-Cummings polariton ladder for 
g/tUQ = 0.2, as well as the evolution of the energy eigenstates of the Rabi model when approaching the strong coupling limit with g ~ coq . 


and quenches on dissipative spin systems. This approach is based on a numerically exact Feynman-Vernon path 
integral approach [52, 53] (see also Sec. 2.4). In Ref. [36], two of the authors together with Z. Ristivojevic and 
R R Orth have extended this approach to dissipative light-matter systems beyond the weak-coupling limit. We 
show extra results in the strong-coupling limit (Fig. 1), in the highly detuned limit A/coq ^ 1, in agreement with 
Ref. [ 140]. The revival probability of the initial state | -F z) ® |0) = |0, -Fz) (with zero photon in the cavity) can be 
computed exactly in the limit A ^ 0, using the shifted b = a^a basis [140]. More precisely, this probability is 
of the form exp—|Z(t)|^ with Z(t) = (g/mo)(e“'“'>'^ — 1) [140], leading to periodic collapses and revivals with 
frequency mo/(27r). The ratio g/cjQ controls the sharpness of these revivals peaks. 

The limit of large spin then referring to the Dicke Hamiltonian (1954) [141] is also solvable through the 
Holstein-Primakoff transformation, exemplifying the possibility of superradiant effects for strong light-matter 
couplings [142]. As mentioned in the introduction, such a superradiant quantum phase transition has been 
recently observed with ultra-cold atoms in an optical cavity, in a non equilibrium setting [5, 6]. The possibility 
of a Dicke quantum spin glass has also been discussed [143]. Some efforts have also questioned the possibility 
of superradiant quantum phase transitions in mesoscopic circuits with Cooper pair boxes [144, 145]. Coupling 
an ensemble of atoms or artificial qubits to light is now commonly realized with various signatures of collective 
effects [22, 23]. We have also shown that entanglement between light and matter in the Dicke model can be 
probed through the properties of the light field [146, 147] or through bipartite fluctuations [148, 149]. 

In the 1960s, Jaynes and Cummings proposed an approximation to the Rabi model [150], referred to as the 
rotating-wave approximation, where the interaction part is simplified into This approximation 

is usually valid when g/too ^ 1 and at small detuning 5 = |a)o — A| -C A. The number of excitations in the 
system (hybrid excitations of light and matter or polaritons), the number of polaritons N = a''a + |(ct 2 + 1)> is 
then conserved as a result of the emergent U(l) s)mimetry (a ^ ae""^ and cr+ ^ We can diagonalize the 

model in the basis |n, -Fz) and \n + 1, —z) where n refers to the number of photons in the system and | ± z) refers 
to the qubit polarization along the z axis. We have the following eigen-energies and eigenstates for iV > 1 (the 
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ground state has energy Eq = 5 /2) : 


Ef,±=NoJo±^^/5^ + Ng\ (2) 

and 

\N+) = a,,\N-l,+z) + l5^\N,-z) (3) 

\N-) = -p^\N-l,+z) + af,\N,-z), 

where = cos(l/2tan“^ _ sin(l/2tan“^ We recover the usual Jajmes-Cummings (JC) ladder 

(see Fig. 1). The anharmonicity of the Jaynes-Cummings ladder produces exotic phenomena such as photon 
blockade, preventing a second photon to enter in the cavity when a first one is also present [151, 152]. The 
Vn spacing of the JC ladder can be used to generate nonclassical states of light. Photon-photon interactions 
occur both in the resonance condition between cavity and qubit and in the dispersive limit, where the light 
frequency and the qubit level separation A are well distinct such that the spin remains in the ground state and 
the dynamics of the photons is described by an effective Hubbard interaction !7(a ' a)^. More precisely, defining 
the small parameter A = g/|a)o — A| and performing an exact unitary transformation on the Hamiltonian results 
in an effective photon-photon interaction U = gA^ [153]. This also produces a shift in the resonator frequency. 
The photon-photon interaction certainly modifies the energy required to add a photon to the cavity depending 
upon the number of photons occupying the cavity. Although this is a higher order effect, the key point is that 
the interaction energy or Kerr energy can be made larger than the cavity line width k and the transport of 
photons through the cavity is blockaded. This photon interaction effect has been observed in various experiments 
(transmission spectrum as a function of frequency or intensity correlation function), for example, in cavity [154] 
and circuit-QED [155, 156, 157, 158]. 

By increasing the light-matter coupling and treating the counter-rotating wave terms perturbatively one 
identifies a Bloch-Siegert shift [159] making the two polariton levels |1—) and |l-l-) repelling each other : when 
entering the strong-coupling limit the level |1—) becomes closer and closer to the ground state (see Fig. 1). In the 
strong-coupling limit, the structure of levels becomes more complicated [56]. In the Rabi lattice situation, features 
reminiscent of the Dicke model appear due to the presence of the counter-rotating wave terms [48]. To get further 
insights, one can project the rabi ground state in the polariton basis [48]. Due to the parity symmetry, only even 
polariton states contribute with different weights, which give the probability distribution of having N polaritons 
in the ground state. This distribution is peaked around some finite value as a result of the counter-rotating terms 
which mix sectors at different polariton numbers, thus acting effectively as a sort of chemical potential. 

In Ref. [36], we used a relatively simple stochastic protocol to address certain non-equilibrium properties and 
quenches on the qubit strongly-coupled to light (see also Sec. 2.4). In particular, we show how to reach an almost 
pure state with one polariton in the branch |1—) by driving the cavity, with an almost zero detuning between light 
and matter. Note that in the limit of large detuning between light and matter, this protocol corresponds to make a 
Ti-pulse on the qubit [160]. Below, we shall discuss further applications. 

2.3. AC driving. Dissipation and Spin-Bath Models 

Let us now discuss dissipation and driving effects in the cavity. First, the cavity line width and driving effects 
from an AC source can be described in a microscopic manner by coupling the cavity to very long transmission 
lines [161, 127] which correspond to a bath of harmonic oscillators [162, 163]. 

Let us envision a configuration where a single cavity is capacitively coupled to a left and a right transmission 
lines (serving as photon baths). The open cavity system coupled to the transmission lines is described through the 
Hamiltonian : 

Hiines= (4) 

j=l,r k k k 

The transmission lines with dispersion = ^1^1 fthe Planck constant h is normalized to unity, for simplicity) are 
then described by the spectral function 

J((u) = 7r^^A^5(a) — cuj.) = 27iaa)exp(—m/o)^), (5) 

k 


5 



and dOj, is an ultraviolet cutoff for photon excitations in the environment (this corresponds to Ohmic dissipation). 
More precisely, defining the symmetric and anti-symmetric combinations of the two transmission lines, only the 
symmetric mode couples to the cavity, and the effective coupling constant a^. turns into 7^^ = aj. Here, 

a^. depends on the properties of the transmission lines, which are supposed to be identical, and ji^r encode the 
capacitive couplings between the cavity and the transmission lines [32, 110]. We relate the bosonic modes in the 
left line bn. with the input signal and the input power P,„ in a usual way, and the output signal is expressed in 
terms of the boson [127]. Using exact equations of motion and the input-output theory [127], we can obtain 
an expression for the transmission and reflection coefficients of a photon in the system [110] : 

r(m) = l-2i^-^J(m)x^^(m) (6) 

Tl +Tr 

t(m) = |t|exp(iy)= 

Tl +Yr 

In this expression, we introduce the exact retarded photon Green’s function XxxM = ~i6(l)([^(0>^(0)])sys 
describing the response of the photon displacement x = (a -I- al)/\/2. The coupling with the environment then 
produces a damping of Ohmic type for the photon propagator [110] : 


Z(<^) 


<^o 

+ (cOqk’ 


(7) 


where k = 2naa> is the effective dissipation induced by a coupling a to the modes of the transmission lines. Notice 
that close to resonance, the dissipation can be described through a complex frequency cjq coQ — iK/l. A damping 
rate linear with frequency is usually typical of long transmission lines which embody a Ohmic environment. Note 
that by treating these dissipation effects in a Hamiltonian manner through a set of bosonic variables already 
constitutes a first many-body effect in these systems. Some dissipation effects can be also described through the 
Lindblad [164] or Bloch-Redfield [165, 166] type equations. Below, we show how many-body physics can be 
engineered from these Ohmic environments. 

Dissipation effects on the spin or the artificial two-level system already results in rich many-body quantum 
impurity models, such as the spin-boson Hamiltonian [52, 53, 167, 168, 169] : 


Hsb = 




A h 1 

+ 2 + 2 2-1 + 


Kl 


( 8 ) 


In the case of Ohmic dissipation (case of Eq. (5)), the physics of this model is very rich showing some exact 
mapping to the anisotropic Rondo model and the classical Ising model with long-range forces [170, 171]; these 
models are now standard models in many-body physics. In particular, this Ohmic spin-boson model yields a 
quantum phase transition for the dissipative parameter ~ 1 by analogy to the ferromagnetic - antiferromagnetic 
quantum phase transition in the Rondo model [172, 173]; here, we use the definition of a given in Eq. (5). It is 
important to mention recent progress in nano circuits [174, 175] towards observing such dissipative quantum 
phase transitions [176, 177, 178, 182, 183]. Rabi oscillations of the qubit also show some incoherent crossover 
for the Toulouse limit [184] corresponding to a dissipative parameter a = 1/2, where the problem becomes 
formally equivalent to a non-interacting resonant level model [185]. Similar quantum phase transitions have 
been predicted for a two-spin system [186, 187, 188] and could be observed in light-matter systems [189]. One 
of the author and B. Coqblin have studied the problem of two-coupled spin-1 impurities [190] and extensions 
to lattice arrays [191]. Non-equilibrium transport at such dissipative quantum phase transitions has also been 
studied theoretically [192] and experimentally [174, 175]. 


2.4. Dissipative Quantum Trajectories, Stochastic Dynamics and Quantum Phase Transition for two spins 

Various methods have been devised to describe quantitatively the dynamics of these spin-boson systems, 
from master equations [193], the non-interacting blip approximation (NIBA) and its extensions [52, 53, 194], 
stochastic wave function approaches [195, 196] and stochastic differential equations [197], time-dependent 
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Numerical Renormalization Group approaches [198, 199, 188], multi-polaron expansion [200], Density Matrix 
Renormalization Group [201] and Matrix Product States [202, 203, 49] for example. These efforts complement 
other methods applied for non-equilibrium dissipative quantum systems based on Monte-Garlo techniques on 
the Keldysh contour [204, 205, 206, 207]. These theoretical developments are very useful owing to important 
progress on the experimental side both in quantum optics, ultra-cold atoms and superconducting qubits, for 
example, by following the real-time dynamics on the Bloch sphere. In the context of superconducting qubits, 
tracking single quantum trajectories has been explored experimentally with a particular focus on the competition 
between continuous weak measurement and driven unitary evolution [16, 17]. Some connections also occur with 
the concept of open quantum walkers [208]. 

At a classical level, the dynamics of dissipative spins in a macroscopic and dissipative environment can be 
described through Bloch equations S(t) = H x S where the dissipation is here included in the noisy random 
magnetic field H = h(cos(^(t),sin<j)(t),0). Classically, the spin then obeys the differential equation S^Ct) = 
—h^ Jjj ds cos((j)(t) — (j)(s))S 2 (s). At a quantum level, as a result of the Ehrenfest theorem, (cj^) obeys a similar 
equation, where we also average over different magnetic field configurations : 


d 

dt 


(crj = -h^ 


f 

Jo 


ds(cos(<j)(t) - (j)(s))CT^(s)). 


(9) 


A usual approximation consists in decoupling the dynamics of the stochastic field and the spin dynamics, neglecting 
completely the back-action of the spin dynamics on the stochastic field evolution. This leads to the well-known 
NIBA approximation [52] (where the noisy phases are then derived from the spin-boson model in Eq. (8) and 
are formally operators [51]). Nevertheless, it is important to mention recent theoretical progress both analytical 
and numerical trying to go beyond these type of approximations [37, 38, 209]. These developments also include 
applications on the driven and dissipative Jaynes-Cummings and Rabi models [36] and networks of spin-boson 
models, leading to many-body stochastic non-equilibrium (mean-field) theories [210, 211]. Lindblad and Bloch- 
Redfield approaches give a complementary point of view [212, 213] and can even include non-Markovian effects 
[197]. Similar approaches have been developed to address heat transfer [194]. 

In Eig. 1, we show one result obtained for the driven Rabi model, by developing a stochastic Schrddinger 
equation approach. The formalism is derived on the path-integral real-time formalism [162] and allows to treat 
Eq. (9) in an exact manner [37, 38]. This formalism has been developed based on previous works in the field 
[39, 40, 214]. In Ref [36], the formalism was developed further in the context of the driven and dissipative Rabi 
model. In Fig. 1, we compute the spin dynamics exactly in the strong-coupling limit, by using (coupled) stochastic 
Gaussian variables to decouple long-range interactions in time induced by the environment (here, the cavity). The 
formalism has also been developed to incorporate dissipation effects, and is also inspired from Ref [214]. In Sec. 
3.1, we show one application in the driven Jaynes-Cummings lattice including dissipation in the cavities. 

In Fig. 2, we show a recent extension of the stochastic Schrddinger equation approach in the context of 
two spins in a ohmic environment [210]. The model is identical to that in Eq. (8), but now one adds an extra 
direct (Ising) coupling Ko'ij.o'zz between spins and the transverse field A acts on each spin separately [188]. 
There are two phases, a paramagnetic phase where the two spins are unpolarized and couple to the bath, and 
a polarized (ferromagnetic) phase where the spins are frozen in a | TT) state, and decouple from the bath : the 
quantum phase transition is of Kosterlitz-Thouless type by analogy to the single-spin problem. The quantum phase 
transition line (characterized by the parameter aj is determined by considering the evolution of the system with 
the two spins being initially in the triplet state |T+) = +z)- For a < a^/2, we observe Rabi oscillations and a 

progressive damping towards the unpolarized equilibrium state. For a^/2 < a< a^, we observe an exponential 
relaxation towards the unpolarized equilibrium state (by analogy to the single-spin problem for a > 1/2). For 
a> a^, the spins stay localized in the initial state |T+) = |+z, + 2 ). This dynamical procedure allows to determine 
quantitatively the position of the dissipative quantum phase transition in the two-spin system. Results for the 
position of the quantum phase transition are also compared with those from the Numerical Renormalization 
Group approach [188] and from Quantum Monte Carlo [215]. 
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A/tj, 


Figure 2 : (Color online) (Left) Phase Diagram of the two spin system obtained from the stochastic Schrddinger (SSE) Equation method for 
two spins [210], in the context of Rabi oscillations after preparing the spin in the initial state |r+) . Here, a is the dissipative parameter of Eq. 
(5), K represents an Ising (direct) coupling between the two spins and A is the transverse field acting on each spin. This model exhibits a 
Kosterlitz-Thouless type transition [187, 188], by analogy to the single-spin problem, and allows to synchronize the spins via the bath as a 
result of the many modes in the environment [188]. The unpolarized phase refers to a paramagnetic phase where the spins are entangled and 
also entangled to the bath, and the polarized phase refers to a phase where the two spins are locked into the (7+) state, which does not evolve 
in time. This quantum phase transition could be realized by coupling the qubits (spins) to a one-dimensional transmission line, Josephson 
junction array, or equivalently to a one-dimensional EEC [226]. (Right) Comparison between the SSE method (blue) [210], the NRG approach 
(red, dashed) [188] and the QMC approach (green, dotted) [215] is shown. There is no direct coupling between the spins [K = 0). 


2.5. Exploration of Josephson-Kondo Circuits with Light 

Josephson physics [216, 217] has been explored with superconductors [218, 219] and quantum circuits 
[220, 129]. Analogous Josephson effects have also been explored with superfluids, for example, in ultra-cold 
atoms [221, 222], in superfluid He [223] and with polariton condensates [224]. 

Here, we review small superconducting circuits in the microwave domain allowing to realize Kondo physics 
with one-dimensional Josephson junction arrays (mimicking long transmission lines) [32, 33] instead of electron 
reservoirs, as in the original Kondo Hamiltonian [26, 28]. The Kondo effect has also been realized in quantum 
dot systems [179, 180, 181]. Excitations in the long transmission lines at equilibrium are zero-point fluctuations 
(bosons instead of electrons in the original Kondo model). Similar geometries have been envisioned in mesoscopic 
rings [182], with Luttinger liquids [183, 32] or one-dimensional BECs (Bose-Einstein condensates) [225, 226]. 
In Fig. 3A, we introduce a typical Josephson-Kondo circuit where micro-wave photons are sent across the system 
[32, 33]. This geometry is also inspired from Refs. [227, 228]. Related circuits with photons have been addressed 
theoretically [34, 35, 229, 230]. Some (fermionic) Kondo type impurity models have also been explored in 
ultra-cold atoms [231, 232]. 

The dissipative parameter in Eq. (5) representing the one-dimensional transmission lines here takes the form 
a ~ 2R/Rq, where R is the resistance of each (superconducting) transmission line and R^ = h/(2e)^ (where 
h = 2 ti within our units) is the quantum of resistance. With the current efforts in nanotechnology one can 
build transmission lines with quite large impedances [174, 175, 233]. In particular, realizing highly dissipative 
transmission lines with Josephson junction arrays is also possible thanks to experimental progress [122]. The 
artificial atom was originally considered to be a double-dot charge qubit at resonance allowing for the coherent 
propagation of Cooper pairs across the system in the micro-wave domain. One could use superconducting qubits 
with long decoherence times [117,118, 119,120, 121,122]. It is also relevant to note recent progress in realizing 
such devices experimentally [234]. 

The Hamiltonian representing the pseudo-spin (representing the double-island charge qubit [235, 236]) 
coupled to the left and right transmission lines takes the form [32] : 

H = + h.c. + ^ 2^'1^1 • (10) 
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The Josephson phases depict the (generalized) superfluid phases in each transmission line in the vicinity of the 
artificial qubit. Note that in the geometry of Fig. 3 (left), only the antisymmetric mode (<f>; — <f>r) couples to 
the transmission line. In the underdamped limit (0 -C ct ^ 0.3) of the spin-boson model introduced earlier, the 
Kondo physics here corresponds to the perfect Rayleigh transmission of a micro-wave photon at the frequency 
corresponding to the Kondo frequency of the system [32]. Such a concept of resonance fluorescence has been 
observed with artificial superconducting qubits [237]. Note that this frequency is distinct from the bare frequency 
of the artificial atom (qubit) as a result of the quantum point-fluctuations in the long transmission lines. More 
precisely, from the two-level perspective, one can define an effective transverse field 


£j{C0S($; -$,)). (11) 

Close to equilibrium, when the input power can be treated in linear response theory, the effective transverse field 
results from known results on the Ohmic spin-boson model [52, 53] : 



The In—corrections which stem from the many-body excitations in the long transmission lines can be understood 
as renormalization effects. This scenario assumes that the high-frequency cutoff co^ of the order of the plasma 
frequency in the transmission lines is much larger than the Josephson frequency. The key point is that the spin 
susceptibility = ~i^^(f){[o'z(f)>o'j,(0)]) takes a form identical to Eq. (6) but with the frequency coq 

being replaced by the Kondo frequency cop and a damping rate y(a)) = ojpJ{cS) which is reminiscent of the 
Korringa-Shiba relation in the spin-boson model [238] and in the Kondo model [28, 239, 240, 241, 243, 242]. 
Note that this damping rate can be measured directly in quantum RC circuits through the concept of quantum of 
resistance [244, 245, 246, 247, 248]. Multichannel effects have also been discussed theoretically [249, 250, 251]. 

To evaluate the elastic transport of a photon through the system, we can still use the input-output theory 
[127] and relate the reflection and transmission amplitudes of a photon [32]. When the averaged input power 
tends to zero, in the underdamped limit, we find that the scattering matrix is unitary and at the Kondo frequency 
the transmission is maximum as a result of J(u)jj)Imx^;^(a)jj) = — 1 (using the conventions of Eq. (6)). The 
phase of the reflected signal also changes from zero to n by varying the asymmetry between the left and right 
couplings with the transmission lines; see Eig. 3. 

Note that inelastic Raman effects nevertheless become important when increasing the drive amplitude or by 
reaching the ultra-strong light-matter coupling (over damped regime), as described in Ref. [33]. Open questions 
concern, for example, the evolution of the Mollow triplet [237] by increasing the drive amplitude and the 
dissipation strength a. Simple arguments allow to conclude the disappearance of the elastic Kondo peak when 
increasing the amplitude of the drive [32]. Similar geometries could be considered to simulate the two-channel 
Kondo effect with light, as in electron systems [252, 253, 254, 255, 256], with finite-frequency measurements 
[257, 258]. 

2.6. Hybrid systems 

Coupling qubits (artificial spins) to a circuit quantum electrodynamics environment allow to explore new 
limits in quantum computation, where entangling distant qubits through light [111] has become possible for 
example [259, 260, 261]. One can also envision the realization of new quantum gates [262] and very robust 
Einstein-Poldolsky-Rosen or Bell’s pairs [263] by analogy to photon systems [264] and Cooper pairs [265]. 
Coupling a real spin-1/2 object to light has also been realized experimentally [266, 267]. Extending the concepts 
and techniques of electrons to quantum optics has also resulted in new developments in quantum opto-electronics 
[268]. 

Recently, experimentalists have accomplished another step further by designing hybrid systems, containing 
electrons with a circuit quantum electrodynamics environment. Eor example, coupling a quantum dot to an 
electromagnetic resonator has been achieved experimentally [112, 113, 114, 115]. Similar systems have been 
built with Nitrogen Vacancy centers [23]. More precisely, the coupling between a quantum dot, a microwave 
resonator and two biased electronic reservoirs was built on chip and signatures of many-body correlations have 


9 



\A/V 


\A^ 




Figure 3 : (LefO Typical Josephson-Kondo circuit introduced in Refs. [32, 33]. In the underdamped limit, where dissipation in the long 
Josephson junction chains is not so important (0 < a C 0.3), one can observe the perfect Rayleigh transmission of a microwave photon 
through the system at the renormalized Kondo frequency toj; given in Eq. (12). This many-body Lamb shift exemplifies renormalization effects. 
Such long transmission lines with tunable resistances can be built with current technology [122]. By adjusting the ratio between the left 
Yl and right couplings between the Josephson junction chains (transmission lines) and the artificial double-dot Cooper-pair box system 
[235, 236], the phase of the reflected signal at the Kondo frequency passes from 0 to rr when changing the ratio Jihr- (Right) Cartoon (note 
the cavity in yellow is much bigger than the mesoscopic (nano-) circuit, as in Fig. 1) : Hybrid system comprising a quantum dot tunnel coupled 
to two (red) electron reservoir leads and capacitively coupled to a one-mode (yellow) cavity via a capacitive coupling, as realized in Ref. 
[112]. In Ref. [110], we have shown that at the renormalized frequency to* (defined in the text), in the non-linear limit of large bias voltages 
applied across the quantum dot, the light is also perfectly transmitted for symmetric couplings y; = yr • 


been observed, for example, in the phase of the microwave transmitted photon as a function of the bias voltage V 
applied across the leads. The system can be modeled by an Anderson-Holstein type Hamiltonian [269, 270] : 

^matter E ^ka^kaa^kaa Vkaiclaa^o + h-C.-J + e E + Unjjii + A[a + a'^J^d'^d^ + coga'a. (13) 

kcra kaa cr a 

Here, and d^ are the annihilation operators of the leads and of the dot respectively with spin u =T, i- In 
addition, = e]^±eV /2 and = dj^d^. The cavity is also coupled to two long transmission lines, similar to Eq. 
(4), that serve to send and receive the micro-wave signals. A similar Hamiltonian has been addressed to study 
transport through molecules with particular modes of vibration [271, 272, 273]. 

One natural question concerns the phase of the transmitted photon signal as a function of the bias voltage 
V across the two electronic reservoirs [110]. Let us consider the quantum limit in the cavity with a few photon 
number. The phase of the transmitted (photon) signal can be computed through the relation tan<p((D) = 
— Re;tf^^(a))/Im/^^(£u) in Eq. (6). By including the cavity-dot coupling, the (light) susceptibility becomes 
modified as ~ tUo/(a)^ — aji^ + (coqK — a)on®(a))). We have treated the light-matter coupling to second 

order in perturbation theory. Interestingly, in the limit of large bias voltage applied across the quantum dot, the 
damping of the light field is only weakly affected by the (fast) motion of the electrons (Im H^CtUo) Re H^Ccuo)), 
then resulting in a phase (p(a)o) which tends to a universal value n/2m the limit of very small dissipation from the 
















cavity environment (tc ^ 0) [110]. Note that this phenomenon occurs both when the quantum dot is at resonance 
and in the Kondo regime. Note that this phenomenon occurs both when the quantum dot is at resonance and in 
the Kondo regime. It is also relevant to note that at the renormalized frequency co* defined through the equation 
co^ — cOq — (UoRen^(a)) = 0, the transmission coefficient exhibits a pronounced peak (see Fig. 2). At large bias 
voltages (D* converges towards coq and the transmission peak approaches unity for symmetric couplings with the 
left and right transmission lines (in Eq. (6)). In addition, the phase of the reflected signal is n by analogy to Fig. 3 
(bottom left) and has been recently experimentally measured in double-quantum-dot geometries in graphene 
settings at low temperatures [274]. It is also relevant to mention other recent efforts from the theoretical point of 
view to quantitatively describe transport in these systems [275], in particular at low-frequency in relation with 
the Korringa-Shiba relation [276]. The strong-coupling limit between light and matter has been also achieved 
experimentally [277]. 

The effect of the electronic non-equilibrium environment on the cavity photon is to induce non linearity, 
damping/friction term and noise term. These terms in general are independent due to current flow and non¬ 
equilibrium effects and can be fully tuned by changing bias and gate voltages. In the limit of large bias, we can 
formally integrate out all the fermionic modes to obtain an explicit expression for the cavity photon effective action 
[110]. From this we can see how these hybrid systems can also serve to simulate novel regimes of the Langevin 
equation and some chaotic properties of the Duffing oscillator [110], which has been realized experimentally [278]. 
Other recent theoretical developments can be found in Refs. [275, 279, 280, 281, 282, 283]. Similar ideas to 
build dissipative arrays have been addressed with hybrid system technology [284]. It is also relevant that coupling 
topological systems such as (one-dimensional) topological superconductors [285, 286, 287] to circuit QED has 
also attracted some attention, both for simulation and measurement point of view [288, 289, 290, 291, 292, 293]. 
AC Josephson effects in these systems have also attracted some attention [294]. Some theoretical efforts have 
also been done to engineer topological Kondo ground states with Majorana quantum boxes [295, 296, 297] and 
multichannel Kondo models [298]. Ideas to implement fermionized photons have also been considered [299]. 
Engineering Majorana fermions with photons has also engendered some attention in the community [300]. 


3. Large Circuits 

3.1. Circuit QED array : Interaction Effects, Effective Chemical Potential and Driving 

In the last few years, some experimental progress were also accomplished towards realizing larger arrays in 
cQED [10, 12]. Experimental efforts have been realized in various geometries [301]. Two cavities can be coupled 
through a capacitive or Josephson coupling [189, 302]. This point has also been carefully addressed in Ref. [124]. 
In particular, a cavity-dimer was realized experimentally [189] allowing to investigate novel non-d 5 mamical steady 
states and entanglement features [303]. 

Assuming that the coupling between cavities (usually in the MHz range) is much smaller than the internal 
frequencies of the cavities, the coupling between cavities turns into an effective hopping term of bosons (photons) 
from one cavity to the neighboring ones. Eor a large quality factor, we may restrict our treatment to a single photon 
mode, similar to the Jaynes-Cummings, Rabi and Dicke Hamiltonians. One model that has been investigated a lot 
theoretically in the literature is the Jaynes-Cummings lattice model, which takes the form [10] : 

H = (ajaj + ajui) - Me// Xi » (14) 

‘ (ci> ‘ 

where describes the Jaynes-Cummings Hamiltonian in each cavity, introduced earlier, and the bracket notation 
(i;J} denotes summation over nearest neighbors (see Eig. 4). One can realize experimentally complex hopping 
tunneling elements J at each junction, by inserting a Josephson junction ring, which allows to break-time reversal 
symmetry and to realize topological phases [61, 62, 64, 65], as described below. Similar quantum nano-circulators 
have also been designed in cQED architectures [63]. This model has attracted some attention, for example, in 
the light of realizing analogues of Mott insulators with polaritons [41, 45, 46, 47]. More precisely, one observes 
two competing effects in these arrays. Eor large values of J, formally the system would tend to form a wave 
delocalized over the full lattice, by analogy to a polariton superfluid [80], whereas for very small values of J, the 
photon blockade in each cavity addressed earlier in Sec. 2.2 will play some role. Theoretical works [41, 46, 47], 


11 


circulator 





t/t^ 

Figure 4 : (Color online) (Left) Example of Superconducting Quantum Electrodynamics Networks, from Refs. [10, 45]. Some experimental 
efforts are currently realized to engineer these artificial photonic arrays operating in the microwave domain [12]. Some other photonic 
(polaritonic) lattices have been built [105, 103]. (Right) Synthetic gauge fields on photon systems can be generated through nano-circulators 
or Josephson junction rings at each Y-junction, on the honeycomb lattice [61]. Efforts are currently developed to realize such nano-circulators 
on-chip [65]. This honeycomb lattice with artificial gauge fields is also dual to a Kagome lattice with particular flux configurations [61, 64], as 
described below. Below, we also show how these gauge fields can be realized through the Floquet theory. (Bottom) We apply the methodology 
described in Ref. [36] and include dissipation effects in each cavity (coq —> tOg ~ tK-/2). The protocol is as follows : we drive the cavities 
to reach a polarition state |1-) in each cavity at a (sufficiently short) time scale t^. Then, we increase the coupling between cavities and 
observe an intermediate state in time, where the standard deviation of the polariton number N = a^a + ^[cr^ -F 1) is negligible (similar to a 
Mott phase). At long times, we identify a dynamical ‘transition’ towards a delocalized state. The parameters are : g/too = 0.02, A/tUg = 0.9, 
kb/tug = 0.1 (Vg is the amplitude of the AC signal on each cavity), k/coq = 10“®, J = Vj(t - tj) (for t > tj) with Vj = 10“^ (in units of 
a>gt“^) and = 311. A similar procedure has been considered to compute the g 2 correlation function in time [304]. 


including ours [45, 51], have focused their attentions on solving the phase diagram at equilibrium in the presence 
of a tunable chemical potential ^. For / = 0, the ground state at small J would correspond to the vacuum in 
each cavity. By driving the system with a quite strong drive amplitude, one can eventually reach for intermediate 
times, a state with one polariton in average in each cavity [36], by analogy to applying a n pulse on the qubits for 
a large detuning between light and matter [160]. By driving cavities individually, an intermediate state analogous 
to a Mott state, was also identified [36]. In Fig. 4, we apply the methodology of Ref. [36], and apply for the 
following protocol: we prepare the cavity in a |1—) polariton state by driving each cavity, then we stop the drive 
and we switch on the coupling between cavities. First, we observe that the effect of dissipation on each cavity 
described by the modification cjq—*cjq — Ik/I (see Eq. (7)) produces some small leakage on the mean polariton 
number per cavity (iV) = (a''a + ^(cr^ +1)), compared to Ref [36]. At longer times, the coupling between cavities 
increases substantially the fluctuations, producing dynamically an extended state for the photons (polaritons). 
The coupling between cavities here is treated at a mean-field level. A possible extension could use the dynamical 
mean-field theory [305]. 

By changing ^ and keeping J small, one could eventually turn the vacuum in each cavity into a polariton 
state in the branch |1—) (by analogy to Fig. 1). Simple energetic arguments predict that this change would 
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occur when Ei_ — = Eq (we use the notations of Sec. 2.2). This result can be made more formal by using a 

mean-field theory and a strong-coupling expansion [41, 45, 46, 47]. In the atomic limit where J is small, one then 
predicts the analogue of Mott-insulating incompressible phases, as observed in ultra-cold atoms [306], where it 
costs a finite energy to change the polariton number [41, 44]. This point has been summarized in various recent 
reviews [10, 11]. By increasing the hopping J, one can build an equivalent of the theory, where ip ~ (u;), 
in order to describe the second-order quantum phase transition between the Mott region of polaritons and the 
superfluid limit [45]. These results, which make a connection with the Bose-Hubbard model [307, 308], have 
been confirmed by numerical methods such as Quantum Monte Carlo [309]. The effect of disorder has also been 
addressed using the Density Matrix Renormalization Group, leading potentially to glassy phases [310]. Disorder 
can also result in Coddington-Chalker type models for light [311]. Weak-localization effects have been observed 
in superconducting circuits [123] .These spin-boson networks can also serve to emulate various quantum spin 
models [45, 43, 312]. It is also relevant to stress the appearance of novel phases in the limit of Rabi lattices, 
beyond the weak Jaynes-Cummings coupling limit; the presence of the counter-rotating wave terms somehow 
acts as an effective chemical potential [48]. The limit of strong light-matter coupling has been achieved with two 
stripline resonators [313]. Dissipation effects in the two resonator problem have been addressed in Ref. [314]. 

Similar to ultra-cold atoms [18, 19], it seems important to be able to tune the effective chemical potential of 
photons to develop further quantum simulation proposals and observe these analogues of Mott phases in the long 
time limit. Recently, an idea to simulate an effective chemical potential for photons was pushed forward in Ref. 
[50], based on a parametric coupling with a bath of the form ?iCos[cOpt)HsB where where 

Bj is a bath operator. Superfluid-Mott transitions using this scheme were also addressed in finite lattices. The Bose- 
Hubbard model was considered which can be achieved for highly detuned cavities [153]. At a more general level, 
considerable attention has been turned towards describing driven cavity arrays in order to realize novel steady 
states. In this context, it is certainly important to develop analytical and numerical tools [315, 316, 317], from 
stochastic non-equilibrium mean-field theories [210, 211] to Matrix Product States for example [49, 318, 319]. 
We have also done some progress to tackle the real-time dynamics of a spin dimer and the quantum Ising model 
in a transverse field coupled to a bosonic environment out of equilibrium [210]. See also recent Refs. [320, 321]. 
Features of the quantum Ising model has been recently observed in cQED [13]. Note that the time-dependent 
numerical renormalization group is also powerful to describe the dissipative two spin system [188]. 

Based on all these theoretical developments, we expect more experimental developments in cQED multi-cavity 
systems [10, 12, 313]. 

3.2. Introducing Topological Phases 

It is also important to emphasize the recent realizations of topological phases with photons in artificial systems 
together with developments in ultra-cold atoms. While in nature this is achieved by a magnetic field coupled to 
charged particles, artificial gauge fields have been realized with ultracold atoms in optical lattices [81, 82] and 
photonic systems [79]. We may think of a topological phase of noninteracting particles as characterized by a bulk 
topological invariant which necessarily implies the existence of edge states protected against backscattering at 
the system boundary. The lure of optical lattices and photonic systems is that band topology, which underlies the 
characterization of quantum Hall-like phases, as well as edge state transport, can be probed. Moreover, a strong 
magnetic field at a suitably chosen filling in the presence of interactions leads to the fractional quantum Hall 
effect. This phenomenon, specific to two-dimensional electron gases, underlies multiple theoretical proposals for 
ultracold atoms [322, 323, 324, 325, 326, 327, 328] or photons [329, 330, 331, 332], but no realization to date. 

Defining (time-dependent) topological invariants in relation with observables is certainly timely for Eloquet 
theories [333, 334]. Measuring topological invariants in quantum fluids can also be done exploiting the Karplus- 
Luttinger anomalous velocity [94, 93, 64]. Another important progress concerns the first realization of a photonic 
topological Eloquet insulator [73]. Rechtsman et al. translated the modulation from the time domain to the spatial 
domain, leading to the experimental demonstration of the photonic analogue of the quantum Hall effect at optical 
frequencies (633nm). Similar developments have been done in micro-cavities where a spin-orbit coupling for 
polaritons has been engineered [76]. Concerning quantum circuits, a topological insulator has been realized 
based on meta-materials, composed of capacitively coupled high-Q inductors [335], and other suggestions have 
been addressed [336]. Other quantum circuit architectures have been proposed [337]. The effect of dissipation in 
driven lattices has also been considered theoretically, in particular, on honeycomb lattice systems [91]. Topological 
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edge states have also been observed on photonic honeycomb lattices [338, 339, 340]. Novel ways of pumping a 
non-equilibrium photonic system were also discussed in Refs. [341, 342]. In fact, these developments are not 
restricted to two dimensions in these quantum fluids by analogy to the topological properties of Bismuth materials 
[69, 70, 71]. 

First, we will review a generic model, the Haldane model, on the honeycomb lattice and its extension to the 
Kagome lattice. Then, we will describe how to implement a Peierls phase for neutral particles based on the Floquet 
mechanism. We will continue with quasi-one-dimensional ladder-type geometries described by a Bose-Hubbard 
type model, related to experiments in ultra-cold atoms and Josephson junction arrays. These systems are also 
related to circuit-QED systems since the Jaynes-Cummings lattice model is described by the same theory as the 
Bose-Hubbard model [45, 51]. 

3.3. Artificial Graphene, Kagome lattice and Haldane Model 

Some experimental works have also succeeded in realizing topological phases through the Haldane model and 
analogous ideas. Here, we review these developments starting from the artificial graphene on the honeycomb 
lattice [105, 343, 344]. The Haldane model takes the form (here, hopping elements are denoted t^ and t 2 by 
analogy to condensed-matter systems) [345] : 

H =-ti^cjcj - t2 ^ +h.c.-F^(-l)'Mcfci. (15) 

{i;j) i 

We fix the effective chemical potential of photons to be zero here. The presence of the term t^ will produce the 
Dirac particles on the graphene lattice and the term t 2 which couples the next nearest neighbor sites corresponds 
to the Haldane term. We also add the Semenoff mass M which has been introduced in high-energy physics [346]; 
for cQED arrays, it corresponds to an adjustment of frequencies in each cavity and in ultra-cold atoms it can 
be engineered through magnetic field gradients, for example. In the continuum limit, we will make an analogy 
with the equation {ij^d^ — M')ip =0 where ip will be a two-component spinor and the Clifford algebra will 
be generated from the Pauli matrices.The configuration of fluxes for the Haldane model is shown in Eig. 5 and 
(pij = (p . The key point here is that the t 2 term opens a gap at the Dirac points by breaking the time-reversal 
symmetry : this produces the emergence of chiral edge modes at a single-particle level. In fact, the system exhibits 
a non-trivial topology as long as M < 3\/3t2 sin cp [345]. 

The precise analogy with the Dirac equation in two spatial dimensions can be seen by introducing the spinor 
ipi = (a^, bi), where a and b correspond to the two sub-lattices [347]. When t 2 = 0 two inequivalent Dirac points 
K and K' naturally emerge in the band structure [348] and are protected by the PT symmetry (Time-reversal and 
inversion S 5 mimetry); for an introduction, see Refs. [347, 349]. The dispersion relation close to the Dirac points in 
that case obeys 4- (Vyk^y, where v,- are the group velocities (below, V; = v), and one can check that 

(PT)H(kXPT)~^ = H(k)*. The implementation of Dirac particles has also been discussed in the cQED context 
[350], in relation with existing experiments. Two linear dispersion bands in the three-dimensional momentum 
space can also intersect at a single degenerate point - the Weyl point. Recently, these three-dimensional Weyl 
fermions have also been identified in photonic crystals [104]. 

The total Hamiltonian can also be written as an effective spin-1/2 particle subject to a magnetic field in the 
wave-vector space : H = L , dkip’’'(k)Hff(k)ip(k) where = —h(k) ■ &, and in the basis of the Pauli matrices 
(a,., <jy, afj the effective magnetic field takes the form : 

h(k) = (ti Z cos(k.aj), tj sin(k.aj), -\-M — 2t2 sin(k.bj)). (16) 

i i i 

Here, and b; are the vectors connecting the nearest-neighbors and the next-nearest neighbors respectively and 
we assume that the phase (p = ti/2 such that the t 2 term in Eq. (15) is purely imaginary. The presence of the t 2 
opens a gap and produces a quantized Chern number similar to the quantum Hall case [345]. In the presence of 
P (T broken), the Berry curvature satisfies JCk) = IF(—k) and the total Berry flux adds up to 27i and the Chern 
number equals one. By analogy with a spin-1/2 particle, the Chern number of the lowest band can also be written 
as : C_ = (l/47r) dkh • (d^h x dy\\), where h is the normalized h vector (here, k lives on the Brillouin zone 
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Figure 5 : Left: Representation of the Haldane model on the honeycomb lattice [345] [top) and spectrum of the system in the topological 
phase [bottom). Right: Kagome lattice with artificial fluxes [64], dual to the honeycomb lattice with Josephson junction rings introduced in 
Ref. [61]. This can be seen as a realization of the Haldane model on the Kagome lattice with no net flux on a parallelogram unit cell. Bottom, 
middle : spectrum of the Kagome lattice for one particular flux configuration (j> = 7i/6 [with the notations of Ref. [64]), showing a flat band 
in the middle of the energy spectrum in addition to chiral edge modes. Right: the local density of states then also shows states trapped on 
hexagonal cells [64], as observed in the experimental setup of Ref. [105]. 


torus whereas h lives on the Bloch sphere). A pedagogical derivation is done in Ref. [349]. An important point is 
that the t 2 contribution is odd under changing k in -k. The bulk-edge correspondence ensures the presence of a 
chiral edge mode in this case similar to the quantum Hall situation [88, 90] (see Fig. 5). A similar quantum Hall 
state without Landau levels can be generated through (long-range) interaction effects [351, 352]. 

In fact, Haldane and Raghu have extended this scheme to the field of gyromagnetic photonic crystals with 
ferrite rods and the magneto-optical Faraday effect then mimics the t 2 term of the Haldane Hamiltonian [87]. 
This theoretical idea has been first realized at MIT in 2009, where a photonic band gap has been observed as 
well as robust chiral edge modes for light in the microwave regime [72]. Note that the experiment has been 
realized with ferrite rods forming a square lattice instead of the honeycomb lattice (see Fig. 6). Photon systems 
are particularly appropriate to produce chiral edge modes, where selecting the right-frequency of the incoming 
wave and observing the local density of states is relatively straightforward [79]. The Haldane Hamiltonian has 
also been achieved in 2014 at ETH Zuerich [353], with fermions in optical lattices, where the t 2 term has been 
simulated using the Floquet point of view and the shaking protocol [77, 78]. Similar experiments have also been 
realized at Hamburg [354]. Preceding experiments on ultra-cold Dirac fermions have managed to move the Dirac 
points by tuning the hopping amplitudes [355, 356]. The Chern number has also been measured in ultra-cold 
atoms [357]. The Zak phase has been experimentally observed [358, 359]. A similar quantum Hall state without 
a uniform magnetic flux was also observed in Bismuth materials [360]. 

It is also important to note that the Haldane model can be extended to more complicated lattices such as the 
Kagome lattice, with three sites (atoms) per unit cell. Following the protocol of Refs. [61, 62, 64], summarized 
in Figs. 4 and 5, one can realize an analogous situation with chiral edge modes of light in cQED. One can map 
the honeycomb lattice with Josephson rings at each link onto a “dual” Kagome lattice with fluxes in the “bow 
ties”. The Kagome model contains particles hopping between nearest-neighbor sites with a complex integral 
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Figure 6 : Top, Left: first experimental observation of chiral edge modes and quantum Hall phases of photons in gyromagnetic photonic 
crystals [72, 79]. Middle, Left: possible Floquet scheme suggested in Ref [381] to realize Peierls phases and quantum Hall phases for photon 
systems without the Faraday effect. Middle Right: Generalization to a quantum Spin Hall state, as observed experimentally in silicon photonics 
at University of Maryland [74, 75]. Top, Right: Spectral function from Ref [395] computed with CDMFT [305, 399] exemplifying that the 
helical edge modes of the Kane-Mele model are robust towards interaction effects. Bottom : Possible phase configuration using the Floquet 
scheme of Eq. (17), producing an (anisotropic) Haldane state with anisotropic couplings (see Eq. (19); here to = n). 


Photons then acquire a phase 3</) around a triangular plaquette and a phase —64) on a honeycomb cell, 
amounting to a zero net flux in a parallelogram unit cell. By adjusting the flux in each Josephson junction ring, 
one can realize a quantum Hall effect without Landau levels [64] (see Fig. 4). This model with three sites per unit 
cell also allows the occurrence of flat bands which corresponds to wave-functions localized on each hexagonal cell 
(see Fig. 5) [64, 361]. The search of flat bands has a long history in condensed-matter physics [362, 363, 364]. 
Note that flat bands have been recently observed at LPN Marcoussis in artificial honeycomb systems [105]. Other 
efforts in polariton systems to engineer and observe flat bands have been done recently both theoretically [365] 
and experimentally [366]. Similar lattices (with fluxes) have been widely investigated theoretically and some 
have been realized in quantum circuits and ultra-cold atoms [367, 368, 369, 370]. The Kagome lattice is also a 
candidate for the exploration of quantum spin liquids (non-ordered Neel states) [371]. Other correlated states 
have been predicted on the Kagome lattice [372, 373, 374, 375, 376, 377, 378, 379]. 

In the presence of disorder, there is however a clear difference between the honeycomb and Kagome lattice. In 
the case of the Kagome lattice, disorder will spread the middle flat band, making it possible for intragap edge 
states to scatter into the bulk. This is consistent with the fact that the Chern number at energies energies close to 
the intragap flat band energy loses its quantization as disorder strength increases [64]. In Fig. 7, we illustrate 
such a difference by computing the Chern number v at a given energy on the disordered Kagome lattice with 
translational symmetry breaking (using non-commutative rules and averaging over disorder configurations) [64]. 
The common characteristic of the flat middle band and the lower band is that as soon as disorder is turned on |v| 
starts decreasing from 1 in the close vicinity of the band edge, as indicated in the plots of Fig. 7. We also notice 
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Figure 7 : Disorder effects on the Kagome lattice. Calculation of the Chern number v in real space at a given energy and at cf) = n/G with 
increasing amounts of phase and on-site disorders. Details of the calculations follow the procedure shown in Ref. [64]. The disorder is a 
white noise (uncorrelated from site to site) disorder. For a disorder amplitude W, a random number between - W12 and W/2 is selected. For 
photonic networks, the left (right) panel depicts the effect of increasing phase (resonator frequency) disorder. Here, to represents the cavity 
frequency. The inset shows the reference bandstructure for the clean translationally invariant system, where Chern numbers corresponding to 
bottom and top bands are -1 and 1, respectively, and the flat middle band has vanishing Chern number. 


that disorder effects have already been controlled and adjusted in superconducting circuits [380]. 

Very recently, another progress was achieved in cQED [100, 101], by observing similar topological effects 
and simulating directly the Haldane-type Hamiltonian H(k^,ky) = v(k^a^ + kyUy) + (M — where 

corresponds to the second-neighbor hopping (tunneling) in a local magnetic field. This Hamiltonian has been 
implemented with one or two superconducting qubits on the Bloch sphere : the BZ on a torus is then simply 
replaced by an integration over the polar and azimuthal angles on the Bloch sphere. A path from north to 
south pole on the Bloch sphere was achieved in time, through microwave pulses. Note that these measurements 
constitute a dynamical (direct) measurement of the Berry phase. These measurements use the general scheme 
proposed in Ref. [102]. Other measurements of the Berry phase were also done in Ref [99]. Motion in a curved 
space will be deflected from a straight trajectory and the state of the system feels a force leading to deviations 

J 7l 

Fg^dO where the 

Berry curvature component Fq^ encodes the information of the force and is related to the spin observable (Uy) 
(the experiments work at constant speed allowing a linear response). The two experiments report a quantized 
Chern number 1 for M < and a topological phase transition for M = m^ where the Chern number drops to 
zero. 

3.4. Topological Quantum Fluids and Spin-Orbit Coupling 

The goal to achieve arbitrary control of photon flows has motivated an increasing recent research on photonic 
crystals, metamaterials and artificial quantum circuits. Here, we first review the theoretical proposal of Ref [381] 
and describe the emergence of protected one-way photon edge states that are robust against disorder, without the 
use of magneto-optical effects [72, 79] (see also Fig. 6, top, left). Using simple protocols, based on particular 
choices of spatial distribution and modulation phases, one can produce an effective magnetic field for photons. 
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leading to a Lorentz force. The Hamiltonian takes the form (see Figure 6, middle left) 

H = + (17) 

i i 

+ V cos(r2t + bj + h.c). 

{i-.j) 

In this scheme, two neighboring resonators have different frequencies. One can simulate various situations by 
adjusting the modulation frequency O. Let us assume that we are close to resonance 0 = 0 )^ — tOg, then we 
can build a simple physical argument based on the rotating-wave approximation. Assuming that V -C O, the 
evolution in time of the photons is controlled by and cug mostly. It is then justified to define the operators 

In the rotating frame, the Hamiltonian is equivalent to charged particles subject to a magnetic field [382, 383] : 

^ + e^'^ucjc;). (18) 

{i;j) 

Using this resonant scheme, then one can produce an equivalent of Peierls phases on a given lattice, for neutral 
particles; here, the nearest-neighbor coupling is dominant. In particular, by adjusting experimentally the phases 
(see Fig. 6, middle left) one can realize for example a uniform magnetic field related to the effective vector 

potential A^f f defined as A^f ^ • dl = . An idea to realize the dynamical coupling in photon (quantum circuit) 

systems has been suggested in Ref [381]. A similar scheme has been used in ultra-cold atoms, based on resonant 
laser-assisted tunneling and tilted lattices [84, 85, 86]. In the context of ultra-cold atoms, tilts can be achieved in 
various ways, from magnetic to gravitational gradients. 

In the case of a uniform net flux, the effective Hamiltonian then exhibits the same properties as a two- 
dimensional gas in the presence of a uniform magnetic field. If the flux is rational, one observes a Hofstadter 
butterfly and chiral edge modes, as achieved in ultra-cold atoms for example [84]. To show this. Fang et al. 
applied a simple Floquet argument [381]. The original Hamiltonian is periodic in time H(t) = H(t -F ^) and the 
eigenstates then take the form where e represents the “quasi-energy” which can be found 

using the equation (idf — H(t))|x(t)) = —e|x(f))- The time-dependent Hamiltonian here takes the particular 
form : H(t) = Hq + -I-where Hq = Hq^ + Hqj, = a),cjc;. The problem can be diagonalized 

using the decomposition | j(t)) = Xt=-oo and using the (a, b) basis, by writing \xn) = \Xn,a) + \Xn,b) 

[384]. Outside the resonant condition, one can apply the Floquet theory, a multi-phase expansion and a 0(l/n) 
expansion [77, 78]. This allows to achieve various flux configurations and spin-orbit couplings [385, 386]. 
In particular, one can realize an effective Haldane-type model, as follows. We expand the effective Floquet 
Hamiltonian in powers of l/H. To first order in l/H, away from the resonance condition, one obtains 

H = ^ ^ sin + 4>ji) b’ib^ + h.c. j . (19) 

We assume that the Hamiltonian also contains a constant nearest-neighbor hopping term (which can be generated 
through a perturbation of the form V cos^(nt -I- <j)ij)). Here, j is the unique nearest neighbor of both sites i and 
k. In Fig. 6 (bottom), we show one application of Floquet theory to realize an anisotropic Haldane model (in 
agreement with the relative signs generated in the terms for the a and b sub-lattices). The generated t 2 terms are 
indeed anisotropic : the t 2 term parallel to the horizontal axis is formally zero. We have numerically checked that 
the Chern number is still quantized and currents propagate at the edges despite the anisotropic choice of phases. 

A generalization of the Haldane model on the honeycomb lattice to spin-1/2 particles has been made by Kane 
and Mele in 2005 [387] (see also Refs. [69, 70, 71]) : 

H = -Yl^lCiaCia+ Xi ^^2VijSlpclcjp+h.C. (20) 

(i;j> {{ Ui))aP 
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The it 2 term (here, purely imaginary) mimics an atomic spin-orbit coupling s^L^. More precisely, on a given 
pair i and j of next-nearest-neighbor sites within the same sub-lattice, we must satisfy = — v^,- = ±1. Similar 
to the Haldane model, the spin-orbit t 2 term opens a gap in the system. The main difference is that now the 
system respects the T-symmetry and is characterized by a different topological invariant. A simple picture of 
the edge states is the Kramers pair (Fig. 6, top right) : a particle (electron) with spin-up propagates to the left 
whereas a particle (electron) with spin-down propagates to the right. This leads to a spin current at the edges 
and to a Z 2 type topological invariant [388, 389] and a spin Chern number [390]. The Z 2 invariant can also be 
expressed in terms of the four time-reversal invariant momenta (points) in the Brillouin zone [391]. Progress 
have been accomplished in ultra-cold atoms, with magnetic field gradients, to observe such a quantum spin Hall 
state [84]. Another important progress from the experimental point of view comes from the discovery of Mercury 
Telluride materials which have allowed to observe a similar physics in two dimensions, leading to the discovery of 
two-dimensional topological insulators [392, 393]. In Refs. [394, 395], together with S. Rachel, W. Wu and W.-M. 
Liu, we have shown that the helical edge modes in the Kane-Mele model are robust towards finite to moderate 
interactions. For free particles, the spectrum of the edge modes can be computed analytically [396]. For recent 
reviews on interaction effects in topological spin-orbit systems and iridate materials, consult for example. Refs. 
[397, 398]. In Fig. 6, we show the spectral functions in the presence of a Hubbard interaction obtained with a 
real-space extension [395] of the Cluster Djmamical Mean-Field Theory (CDMFT) [305, 399]. We have reached a 
similar conclusion using the Slave-Rotor approach [394, 400]. It is also interesting to note the possibility of chiral 
bogoliubons in a EEC system [401]. 

Recently, a spinful version has been realized in photon systems by Hafezi et al. [74, 75] based on a Silicon 
architecture (see Fig. 6, middle) : the two spin-polarizations correspond to the two possible circulations (clockwise 
or anti-clockwise) in a given resonator. The phase between resonators are carefully adjusted to follow the scheme 
of Fig. 6. The main difference with the usual quantum spin Hall effect is that here the two spin polarizations 
couple when T is not broken. To decouple the two spin polarizations, one requires to break the T symmetry; in 
this sense, one realizes two independent copies of a Harper-Hofstadter model; as a result, topological protection 
is only partial and limited to devices with no back-scattering. The measurement of a topological invariant in 
these systems has been suggested [91, 92]. An optical spin-Hall effect in photonic graphene and other spin-orbit 
topological phenoma have also been suggested [402, 403]. 

3.5. Josephson Effect and Chiral Bosonic Phases 

Here, we discuss interaction effects on topological phases. Then, we give a specific example related to bosons 
systems which can be realized in Josephson junctions, ultra-cold atoms or photon systems. 

The effect of interactions on topological phases have been addressed in various contexts [404, 405, 406]. 
The stability of topological insulators, for example, towards disorder has been thoroughly analyzed from a 
theoretical point of view [407, 408, 64]. For problems such as the quantum spin Hall state, one can apply various 
analytical arguments to answer these types of questions [387, 394]. The stability of topological insulators and 
superconductors is in general ensured perturbatively due to the gap protection, for fermion systems. Extending the 
notion of topological invariants to interacting fermion systems has also attracted some attention in the community 
in relation with numerical developments [409, 410, 411, 412]. 

Eor boson systems, the situation is different since the bosons can condense in a Bose-Einstein tjqje manner and 
occupy macroscopically the bottom of the lowest Bloch band. Defining topology in these many-body boson systems 
has attracted attention recently. Here, we review two examples of non-trivial effects in chiral bosonic systems, 
in relation with current experiments [353, 413]. The breaking of time-reversal symmetry via the spontaneous 
formation of chiral order has been observed in Ref. [414]. The physics here is associated with Meissner physics. A 
connection to quantum Hall phases in ladder systems will be addressed. 

The first effect concerns a situation related to Eqs. (17) and (18), which is obtained from the time-dependent 
Hamiltonian at the resonance condition. Here we follow the recent experiment realized at Muenich on a ladder- 
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Figure 8 : DMRG and Exact Diagonalization results on small ladders with bosons, with 8 and 16 rungs, showing that the Meissner currents 
subsist in a Mott phase with total density 1, as predicted by the theory [419, 420]. Here, the parameters are fixed to : formally U —> +oo 
(hard-core bosons), V± = 4t, t = 1, t^ = 1 and the flux is fixed to ;^ = 7t/25 using the notations of Ref. [420]. Other careful DMRG analyses 
of this model have been performed in Refs. [423, 424, 425]. 


type geometry [413], and focus on the effective (time-independent) Hamiltonian : 

H = -t ^ (21) 

ia i ia 

a= 1,2 denotes the two chain index, I is the lattice spacing in both directions (parallel and perpendicular to 
the chains) and we have considered a general gauge field acting on the boson operators of each chain. We have 
also added a chemical potential which is justified in ultra-cold atoms, Josephson junction arrays and can be 
simulated in photon networks, as discussed in Sec. 3.1. Here, we also consider the effect of on-site repulsion and 
repulsion between chains, nj„ = Uj'^Uia- The ty term represents a Josephson coupling between chains favoring 
(superfluid) phase-locking between the chains and the occurrence of a superfluid. The artificial flux phases 
can be explicitly derived from a Floquet approach [415]. The application of a small (synthetic) magnetic field 
perpendicular to the ladder plane then will result in the Meissner effect as in superconductors. We can describe 
this Meissner effect as follows [416, 417]. Let us introduce the explicit bosonic representation a,a = -/ne'®'”, 
where n represents the mean density in each chain, describing the (quasi-) superfluid in each chain. The Josephson 
term can be re-written as —t^ncosQA^i + ®li ^2i)- The current (density) operators at point i take the form : 
_)'ll = it(— -I- + h.c. and jj_ = —+ h.c.. The key point is that for strong values 

of ty, the superfluid phases are pinned such that ®li ^2i -H lAj_i = 0 and therefore we obtain a Meissner effect, 
meaning that (;[l) = 0 in the bulk (see Fig. 6). Furthermore, in the small field limit, we may expand to get 

0||) = (Ji) - ih) = -2tn ifij. (22) 

Here, we have defined the phase around a plaquette, = (A^^ + (An Then, this corresponds 

to the usual Meissner current, here defined on a lattice : a chiral current moves in opposite direction than the 
applied external flux. This is an example of chiral superfluid phase in low dimensions which has been observed 
experimentally [413]. This argument is in fact general since the tj_ hopping term is a relevant perturbation under 
Renormalization Group arguments in the superfluid phase; one can also apply bosonization techniques explicitly 
[418]. By increasing the magnetic flux, one sees the emergence of vortices; similarly to type II superconductors. 
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there is a critical field strength. Another nice phenomenon is that these Meissner currents can subsist when 
entering into a Mott phase. For example, for one boson per rung in average, as shown in Refs. [419, 420], by 
adjusting the parameters U and V]^, one can enter a Mott state where the total density in each rung is fixed but 
still Meissner currents develop : this is an example of chiral Mott phase which breaks time-reversal symmetry. 
In Fig. 8, we show that the Meissner currents in the Mott phase can be already observed for relatively small 
systems (with 8 or 16 rungs for example; these results are obtained with the Density matrix Renormalization 
Group (DMRG) and Exact Diagonalization). More results have been confirmed independently based on the DMRG 
[423, 424]. Other chiral phases have been predicted for other filling factors [425, 426, 427]. Note that Laughlin 
type phases in ladder systems [428, 429, 430] characterized by chiral edge states are also possible by adjusting 
densities and the magnetic flux, as shown in Ref. [420]. These correlated phases can also be realized in quantum 
circuits such as Josephson junction ladders [420]. Here, it is also important to emphasize other experimental 
effort in ultra-cold atoms, with boson and fermion systems [421, 422]. 

Another example of chiral bosonic phases has been be identified through the Haldane model on the honeycomb 
lattice discussed in Sec. 3.2, with bosons. Extensive results have been obtained through a combination of bosonic 
DMFT, Exact Diagonalization and analytical arguments [431]. Here, we review simple physical results. One 
can easily identify the presence of two superfluid phases in the phase diagram. For t 2 = 0 (here, we follow the 
notations of Sec. 3.2) all the bosons condense at the center of the Brillouin zone, at zero momentum. On the other 
hand, for t^ = 0, the model turns into two decoupled triangular lattices, and the bosons condense at the two 
inequivalent corners of the Brillouin zone K and K'. Whereas the superfluid phase at t 2 = 0 will exhibit Meissner 
physics in response to the Haldane phases = 7i/2 between the next-nearest-neighbor bonds, the phase at large 
t 2 constitutes a novel chiral superfluid phase. More precisely, as an order parameter, one can define the average 
current (j^) coupling a pair of next-nearest-neighbors which belong to the same sub-lattice. In the superfluid phase 
with t 2 = 0, the Meissner response to the phase = —cpij = —nil is {j,j) = —2nt2lm(exp(—i7r/2)) = 2nt2. 
In the chiral superfiuid phase, one rather observes a sign change of the response (j^) = —nt 2 . This implies that 
the quantum phase transition between the two phases is a first-order-phase transition. A thorough analysis has 
been performed to build the complete phase diagram [431]. The properties of the Mott phase, emerging for one 
particle per site, are also interesting; the Mott phase exhibits at high energy topological quasiparticle and quasi 
hole excitations visible in the spectral function. Similar chiral superfluid phases have been predicted on the square 
lattice [432]. 

We note that the emergence of chiral currents have been predicted in other Hamiltonians, for example, in 
doped Mott insulators [433, 434], high-Tc superconductors [435] and other ladder models [436, 437, 438]. 

3.6. Quantum Link Model in superconducting circuits 

Substantial theoretical effort has gone into implementing gauge theories in ultra cold atoms and quantum 
circuits [439]. Superconducting circuits have also attracted some attention in the context of topological quantum 
computing [440] and Quantum Error Correction codes [441], with some experimental efforts in superconducting 
circuits [442]. There are also current theoretical efforts to suggest Majorana surface codes with superconducting 
wires (circuits) [443, 444]. One can provide a useful scheme to think about these gauge fields with potential 
realizations in superconducting circuits, using loops of Josephson junctions as described in Ref. [445]. One can 
build a Quantum link model. Let us consider a square lattice and two nearest neighboring sites connected with the 
phase ifiij = Jl dl • A. The flux variable is then defined as L/^ = exp(i(py) e 1/(1), and the electric field (conjugate) 
variable connecting this link is then defined as = —id/dpij. These variables are operators satisfying the 
angular momentum type algebra [£,^,1/;^] = Uij and [E^^, I/l] = —LfJj', Uij and UL then act as lowering and 
raising operators of the electric field e^, such that 

Eij\eij) = eij\eij). (23) 

Here, eigenvalues of the electric field are integer or half-integer. One can write down the Hamiltonian, by 
introducing the flux per plaquette. For example, for a square type geometry, the flux per plaquette reads 
I/q = = exp(i4>) where the magnetic flux per plaquette is 4> = -F -F The dynamics 
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of the gauge fields here is described by the Hamiltonian [445] : 

« = + C24) 

<i;j) □ 

This Hamiltonian corresponds to a Maxwell-type Hamiltonian on a lattice (the second term in cos <f> is related 
to Josephson physics). Lattice gauge transformations keeping the Hamiltonian invariant have been discussed 
explicitly in Ref. [445]. One can envision different Quantum link Models. A class of models making links with 
spin-1/2 Hamiltonians and dimer models is defined as follows. If = ±1/2, then one can define spin-up and 
spin-down states through (£; = Sf,, I/;,- = S/) and (£;; = —= S/[), where Ef- = = 1/4. In this 

language, the Uq represent ring-exchange terms In a neutral sub-space of the Hilbert space we 

have the same number of spin-up and spin-down states. Then, one can also engineer a Rokhsar-Kivelson type 
interaction between dimers [446] 


H = -J 2 ([/□ ± 1/^ - KUn + U'af) • (25) 

□ 

This model has attracted some attention in condensed-matter physics, in relation with confined-deconfined phase 
transitions. Theoretical questions about gauge theories are also relevant [447]. Generalizations to non-Abelian 
situations might lead to more complex gauge theories. Some applications of quantum link models to quantum 
chromodynamics have been suggested [448]. Other gauge theories may have applications for quantum computing 
and the toric code [449]. In fact, theorists have recently shown that linear circuits can already mimic non-trivial 
hopping matrix elements, in relation with current technology, allowing potentially to realize non-Abelian gauge 
fields. In particular, loops of “triangle-like capacitors” have been considered as a potential simulator of a spin- 
doubled Azbel-Hofstadter [383, 450] model, with 1/3 magnetic flux per plaquette [337]. This model is related to 
the time-reversal invariant Hofstadter model [451, 452, 453, 454]. Non-abelian models have also been considered 
in ladder systems [455]. 

4. Conclusion 

We have summarized several recent developments, both theoretical and experimental, to realize many-body 
physics with photons in superconducting Quantum Electrodynamics (QED) Networks. At the theoretical level, 
we have outlined analytical and numerical approaches to describe quantitatively these open, dissipative and 
many-body quantum systems beyond the equilibrium limit. We hope that this review will also serve to inspire 
novel experimental developments, related for example to the realization of Rondo physics with light in Josephson 
junction arrays [32, 33] or to the realization of the Jaynes-Cummings lattice [10, 41, 43, 44, 45, 46, 47] or 
Rabi model [48, 49] with superconducting circuits. Some experimental efforts already exist for a few cavities 
[189, 301, 302]. At a more general level, one can also anticipate closer links with propagating light in bulk 
nonlinear media [457] and with phononics [458]. Dissipation has also been shown to be efficient to induce 
topological states of matter [459]. 

We have presented results on the driven Jaynes-Cummings lattice from a stochastic Schrodinger equation 
approach [36, 37, 38, 39, 40] including dissipation effects in the cavities. The coupling between cavities is treated 
at a mean-field level. A more rigorous dynamical-mean field theory approach could be, in principle, built [305] 
by analogy with recent results on dissipative spin chains with long-range forces [210]. We have also extended 
stochastic approaches for two-spin dimer problems [210], which have attracted some attention both theoretically 
[303] and experimentally [189]. On the experimental side, hybrid mesoscopic circuits [111, 112, 113, 114, 115] 
comprising quantum dots and photons will hopefully confirm our theoretical predictions [110]. Regarding the 
implementation of synthetic gauge fields and spin-orbit couplings in photon systems, experimentalists have 
already managed to observe topological phases of light [72, 73, 74, 75, 76, 79], in synchronization with the 
developments in ultra-cold systems [81, 82, 83, 84, 85, 86]. Using the Floquet approach, we have shown that the 
Haldane model on the honeycomb lattice can be engineered from a perturbation periodic in time. The Haldane 
model has been recently implemented with ultra-cold fermions [353]. A similar scheme has been considered 
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to simulate a Kane-Mele model for photons [74, 75]. Josephson junction loops can serve to realize artificial 
gauge fields and break time-reversal symmetry [61, 62, 64, 63, 65]. Applying the Density Matrix Renormalization 
Group and Exact Diagonalization (on small systems), we have checked that Meissner currents can develop in a 
Mott insulating phase, in a ladder geometry with synthetic gauge fields and Josephson physics [419, 420, 427]. 
Such ladder geometries have been built with ultra-cold bosonic atoms [413] based on the theoretical work 
of Ref. [416]. Experimental accomplishments with lumped elements [335] and small superconducting circuits 
[100, 101] pave the way towards observing topological phases in superconducting Quantum Electrodynamics 
(QED) networks. Entanglement properties in such quantum Hall phases could be probed through charge noise 
[148, 149, 460, 461,462]. 
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